First Principles Study on the origin of large 

thermopower in hole doped LaRh03 and CuRh02 
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' (\ ' Abstract. Based on first principles calculations, we study the origin of the large 

&> , thermopower in Ni-doped LaRhOa and Mg-doped CuRh02 We calculate the band 

structure and construct the maximally localized Wannier functions, from which a 
tightbinding Hamiltonian is obtained. The Seebeck coefficient is calculated within the 
Boltzmann's equation approach using this effective Hamiltonian. For LaRh03, we find 
that the Seebeck coefficient remains nearly constant within a large hole concentration 
range, which is consistent with the experimental observation. For CuRh02, the overall 
temperature dependence of the calculated Seebeck coefficient is in excellent agreement 
with the experiment. The origin of the large thermopower is discussed. 
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1. Introduction 

The discovery of large thermopower in Na x Co02 flj and the findings in 
cobaltates/cobaltites[H [7J [8] and rhodates[9|, [10] that followed have brought up 
an interesting possibility of finding good thermoelectric materials that have relatively 
low resistivity. We have recently proposed that the "pudding-mold" type band is the ori- 
gin for the coexistence of the large thermopower and the low resistivity in this material. 
[2]. Let us first summarize our idea. Using the Boltzmann's equation, the thermopower 
is given as 

S = ^K^Ki (1) 

where e(< 0) is the electron charge, T is the temperature, tensors Ko and Ki are given 
by 

dm 



k 



de 



-(*) 



(e(fc)-Ai)". (2) 



Here, e(fe) is the band dispersion, v(k) = V^e(fc) is the group velocity, r(k) is the 
quasiparticle lifetime, f(e) is the Fermi distribution function, and ft is the chemical 
potential. Hereafter, we simply refer to (K n ) xx as K n , and S xx = (l/eT)(Ki/K ) (for 
diagonal K ) as S. Using K , conductivity can be given as a xx = e 2 K = a — 1/p. 
Roughly speaking for a constant r, K ~ E'(t>^ + v B ),K\ ~ {k,BT)Y/{v B — v A ) (apart 
from a constant factor) stand, where E' is a summation over the states in the range of 
\s(k) — (x\ <~ ksT, and va and vb are typical velocities for the states above and below 
/i, respectively If we consider a band that has a somewhat flat portion at the top (or 
the bottom), which sharply bends into a highly dispersive portion below (above). We 
will refer to this band structure as the "pudding mold" type. For this type of band with 
\i sitting near the bending point, v\ ^> v B holds for high enough temperature, so that 
the cancellation in K x is less effective, resulting in \K\\ ~ {k B T)T!v\ and K ~ E'it^, 
and thus large |>S| ~ 0{k B /\e\) ~ O(100)/xV/K. Moreover, the large va and the large 
FS results in a large K oc a as well, being able to give a large power factor S 2 /p, which 
is important for device applications. 

In the present study, we focus on a possibly related rhodate LaRh0 3 j3] with Ni 
doping and CURI1O2 [T5| fT6] with Mg doping, in which large thermopower has been 
observed. In Ni doped LaRh0 3 , the Seebeck coefficient at 300K steeply decreases up 
to the Ni content of x — 0.05, but then stays around 100/xV/K up to about x = 0.3. 
On the other hand, the conductivity monotonically grows, resulting in a monotonically 
increasing power factor (see Figj2]). For CuRh02, the two existing experiments give 
different results. In ref. [16], the Seebeck coefficient is found to be nearly independent of 
doping, while it decreases with doping in ref. [T5] . 
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Figure 1. The band structure of the distorted (a) and the ideal (b) structure the 
tightbinding mode (solid lines) together with the LDA band calculation results (dotted) 
are shown. 



2. Method 



LaRh03 has an orthorhombic structure, which is distorted to some extent from the 
ideal cubic perovskite structure. The experimentally determined lattice constants are 
a = 5.5242(12), b = 5.7005(12) and c = 7.8968(17) A HI]. For comparison, we also 
calculate the band structure for the ideal cubic perovskite structure, where the lattice 
constant is taken as a = 3.940A.[14j CuRh02 has an delafossite structure whose 
experimental lattice constants are a = 5.810910, c = 32.437162A. We have obtained 
the band structure of these materials with the Quantum-ESPRESSO package [T2]. We 
then construct the maximally localized Wannier functions (MLWFs)[13j for the energy 
window — 1.75eV < e k — E F < — 0.64el^ for the ideal structure of LaRh0 3 , — 1.8eV < 
€k — E F < 0.5e\^ for the distorted structure of LaRh0 3 and — lOeV < e^ — E F < 4eV for 
CuRh02, where e& is the eigenenergy of the Bloch states and Ep the Fermi energy. With 
these effective hoppings and on-site energies, the tight-binding Hamiltonian is obtained, 
and finally the Seebeck coefficient is calculated using eq.flTJ. The doping concentration 
x is assumed to be equal to the hole concentration, and a rigid band is assumed. 



3. Results and Discussions 



We first present results for LaRh0 3 . The calculated band structure for the distorted 
structure of LaRh0 3 is shown in Fig{T] along with that for the ideal structure. The 
tightbinding model Hamiltonian for the distorted structure consists of 12 bands (4 
Rh per unit cell), while the model for the ideal structure contains three ti g bands. 
The calculated Seebeck coefficient at 300K is shown in FigfS] as a function of hole 
concentration together with the experimental result. [3] Here we assume that the hole 
concentration n^ is equal to the Ni content. It can be seen that the Seebeck coefficient 
steeply decreases with doping with n h < 0.05, but stays nearly constant for n h > 0.1 
for the distorted structure in particular. As a result the (normalized) power factor 
monotonically grows with doping, which is at least in qualitative agreement with the 
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Figure 2. (a) Seebeck coefficient for the distorted(green) and the ideal(blue) structure. 
The red line is the experimental result. [3] (b) Power factor of the distorted structure 
normalized at n^ — 0. The red line is the experimental data. [3] 
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Figure 3. DOS of distorted (green) and ideal(blue) structures. The right panel is a 
blow up of the left. 
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Figure 4. Tightbinding band obtained via maximally localized Wannier orbitals 
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Figure 5. Calculation result of the Seebeck coefficient for CuRhC>2 with the hole 
concentration of n^ = 0.05 and n^ = 0.1. The experimental data are from ref.[15j 
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experimental observation. 

Now, in order to understand this peculiar hole concentration dependence of the 
Seebeck coefficient, we now turn to the density of states (DOS). The comparison of the 
DOS between the two structures is shown in Figj3l The DOS at the band top is larger 
for the ideal case since the three bands are degenerate. Thus, for low doping, Ep stays 
closer to the band top for the ideal structure, resulting in a larger Seebeck coefficient. 
This is a typical example where the multiplicity of the bands lead to an enhanced 
thermopower, i.e., the larger the number of bands, the closer the Ep to the band top. 
In the case of the distorted structure, as the hole concentration increases, Ep lowers 
and hits the portion of the band with a large DOS (Figj3left). Therefore, Ep hardly 
moves with doping, resulting in a slow decrease of the Seebeck coefficient. A large DOS 
region lies in a lower energy regime in the ideal structure, and therefore the Seebeck 
coefficient continues to decrease with doping (up to a larger doping concentration) . The 
large Seebeck coefficient of about 80/iV/K in the distorted structure can be considered 
as due to the flatness of the top of the bands (around the T point), i.e., the pudding 
mold type band. 

We now move on to CuRh02- The calculated band structure is shown in FigJU 
Around the T point, there is again a pudding mold type band, whose top is very flat. 
The calculation result of the Seebeck coefficient as a function of temperature is shown in 
Figj5] We find excellent agreement with the experiment in ref. [15] in a wide temperature 
range and for the Mg content x = 0.05 and x = 0.1. On the other hand, in ref.jl6], the 
Seebeck coefficient is nearly independent on x. The origin of the discrepancy between 
this experiment and the present calculation remain as a future problem. 

4. Conclusion 

To conclude, we have studied the origin of the large thermopower in LaRhOs and 
CuRh02- From the first principles band calculation results, a tightbinding model is 
obtained via the maximally localized Wannier orbitals, and the Seebeck coefficient is 
calculated using the tightbinding model. In both materials, the large value of the 
Seebeck coefficient can be considered as due to the flatness of the top of the bands i.e., 
the pudding mold type band. For LaRhOa in particular, the Seebeck coefficient barely 
decreases for the hole concentration of n^ > 0.1 in agreement with the experiment, 
which we attribute to the peculiar uprise of the DOS near the band top. 

The authors acknowledge I. Terasaki, S. Shibasaki, and M. Nohara for valuable 
discussions and providing us the experimental data. Numerical calculation has been 
done at the Supercomputing Center, ISSP, University of Tokyo. This study has been 
supported by the Grants in Aid for Scientific Research from the MEXT of Japan and 
from JSPS. 
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